sort 序号 年份
sort 年份 序号

tab 年份, gen(年份)

tab 序号, gen(序号)

global year "年份2007-年份2021" 
 
gen ln二手房=log(二手房)
gen ln二手占比=log(二手占比)
gen ln收入=log(收入)
gen ln赤字=log(赤字)
gen ln城市投资=log(城市投资)
gen ln城投债=log(城投债)
gen ln住宅面积=log(住宅面积)
gen ln销售面积=log(销售面积)
gen ln住宅投资=log(住宅投资)
gen ln医院=log(医院)
gen ln人口密度=log(人口密度)
gen lngdp=log(gdp)
gen ln垃圾=log(垃圾)
gen ln绿化=log(绿化)
gen ln财政收入=log(财政收入)
gen ln财政支出=log(财政支出)
gen ln中小学=log(中小学)

sum 

  reg ln二手房 ln财政收入 ln收入 ln赤字 ln城市投资,r 
  est store OLS
   predict e, res
   sum e, de  

//模型1,OLS

*-Test sig_w = sig_u = 0   
*- MLE of OLS model
constraint define 1 [sigma_u]_cons = 0
constraint define 2 [sigma_w]_cons = 0  
local xx " ln财政收入  ln赤字 ln收入 ln城市投资"
ml model lf SFA2tier_KP05_extend_ll (ln二手房: ln二手房 = `xx') ///
          (sigma_v:) (sigma_u:) (sigma_w:), constraint(1 2)
qui ml check
qui ml search
ml max
est store m0     //模型2



 *--------------  
 *-two-tier SFA
   
    
   SFA2tier ln二手房 ln财政收入 ln收入 ln赤字  ln城市投资 ln住宅面积 ln住宅投资 ln医院 ln人口密度 lngdp ln垃圾 ln绿化, search     
   version 10.0
   ml max, difficult 
   est store m1     //模型3
   
   
    SFA2tier ln二手房 ln财政收入 ln收入 ln赤字 ln城市投资 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln垃圾 ln绿化  年份2007-年份2021            ///
            , search   
   version 10.0
   ml max, difficult 
   est store m3  //模型4
  
	
   
   SFA2tier ln二手房 ln财政收入 ln收入 ln赤字  ln城市投资 ln住宅面积 ln住宅投资 ln医院 ln人口密度 lngdp ln垃圾 ln绿化  ///
            序号1 序号2 序号3 序号4 序号5 序号6 序号7 序号8 序号9 序号10 序号11 序号12 序号13 序号14 序号15 ///
			序号16 序号17 序号18 序号19 序号20 序号21 序号22 序号23 序号24 序号25 序号26 序号27 ///
			序号28  序号29  序号30 ///
            , search   
   version 10.0
   ml max, difficult 
   est store m2    //模型5
   


   
   SFA2tier ln二手房 ln财政收入 ln收入 ln赤字 ln城市投资 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化 ln垃圾  年份2007-年份2021            ///
			序号1  序号2 序号3  序号4  序号5 序号6  序号7  序号8 序号9 序号10 序号11 序号12 序号14  序号17 ///
			 序号21 序号22  序号23  序号24  序号27  序号28 序号30     ///
            , search   
   version 10.0
   ml max, difficult 
   est store m4    //模型6
  
   SFA2tier_sigs 
   
   
*------------                 
*-结果汇总
*------------ 

 *-屏幕呈现
   local mm "OLS m0 m1 m2 m3 m4 "  
   esttab `mm', mtitle(`mm') b(%6.3f) compress  ///
           nogaps star(* 0.1 ** 0.05 *** 0.01)  ///
           s(N r2_a ll) sfmt(%8.2f)
		   
////////////////////////////////////

local model m1 m2 m3 m4
local k : word count `model'
mat LR2 = J(`=(`k'-1)*2',`k',.)
local mm "`model'"
forvalues i=1/`k'{
   gettoken v mm: mm
   local ii = `i'*2-1
   local j = `i'+1
   foreach subm of local mm{
     qui lrtest `v' `subm'
     mat LR2[`ii',`j++'] = (r(chi2) \ r(p))
   }
}
mat colnames LR2 = `model'
mat rownames LR2 = chi2 p-value chi2 p-value chi2 p-value 

//////////////////////////////////
local model m1 m2 m3 m4 
local k : word count `model'
mat LR1 = J(2,`k',0)
local j = 1
foreach md of local model{
   qui lrtest m1 `md'
   mat LR1[1,`j++'] = (r(chi2) \ r(p))
}
mat colnames LR1 = `model'
mat rownames LR1 = chi2 p-value
dis in g "LR test based on Model(0)"
logout, save(Result\Table02_LRa) excel replace:   ///
    mat list LR2, format(%6.3f) 

 *---------------  
*-估计单边效应 
*---------------  
  drop *_hat* *_diff* pat doc doc_pat
  SFA2tier_eff
  replace uw_diff_exp = -uw_diff_exp 
  
  foreach v of varlist u_hat w_hat uw_diff *_exp{
     replace `v' = `v'*100
  }
 
   *-样本总体
     tabstat u_hat_exp w_hat_exp uw_diff_exp,     ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s) save
    tabstatmat all
    mat all = all'
    mat list all, title(full sample)
	
    
   *-年度统计
      bysort 年份: tabstat u_hat_exp w_hat_exp uw_diff_exp, ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s) 
       tabstat  uw_diff_exp, by(年份) ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s)    
  
                        
   *-urban: 东中西
     bysort 地区: tabstat u_hat_exp w_hat_exp uw_diff_exp, /// 
             s(mean sd p25 p50 p75) f(%6.2f) c(s)   

			 
			 
 *-直方图  u的图
   set scheme s1mono
   histogram u_hat_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("供给者 (%)")  ///
       xscale(titlegap(3) outergap(-2))   
       graph export Result\u_his.wmf, replace

   
   histogram w_hat_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("需求者 (%)")  ///
       xscale(titlegap(3) outergap(-2)) 
       graph export Result\w_his.wmf, replace   /// w的图
	   
	   
	   * replace uw_diff_exp = -uw_diff_exp
   histogram uw_diff_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("净剩余 (%)")   ///
       xscale(titlegap(3) outergap(-2))    
       graph export Result\uw_his.wmf, replace  /// 图3
	   
	   
	   
			 
 *稳健性
 
 
 
   cap program drop SFA2tier    
   SFA2tier ln税基 ln财政收入 ln收入 ln赤字  ln城市投资 ln住宅面积 ln住宅投资 ln医院 ln人口密度 lngdp ln垃圾 ln绿化, search     
   version 10.0
   ml max, difficult 
   
    SFA2tier ln税基 ln财政收入 ln收入 ln赤字 ln城市投资 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化 ln垃圾  年份2007-年份2021           ///
			序号1  序号2 序号3  序号4  序号5 序号6  序号7  序号8 序号9 序号10  序号11  序号13 序号14  序号12  ///
			序号16                    ///
            , search   
   version 10.0
   ml max, difficult 
   
   SFA2tier ln二手占比 ln财政收入 ln收入 ln赤字 ln城市投资 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化 ln垃圾  年份2017-年份2021            ///
			序号2 序号4  序号6 序号7  序号8 序号10  序号11  序号12  序号14  序号16 序号20   序号21    ///
			序号23 序号25  序号27 序号1  序号3 序号17       ///
            , search   
   version 10.0
   ml max, difficult  
   
  SFA2tier ln二手房 ln财政收入 ln收入 ln赤字 ln城投债 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化 ln垃圾  年份2016-年份2021           ///
			序号4 序号5 序号7 序号10  序号11 序号12  序号13 序号14  序号16 序号20  序号21 序号27    ///
            , search   
   version 10.0
   ml max, difficult  
   
  SFA2tier ln税基 ln财政收入 ln收入 ln赤字 ln城投债 ln住宅面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化 ln垃圾  年份2016-年份2021            ///
			序号1 序号2  序号7 序号9  序号10 序号12  序号28   ///
            , search   
   version 10.0
   ml max, difficult  

 
	   
